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ABSTRACT 

This work is a continuation of plasma oscillation studies 
made in rare gas plasmas extending over the last four years at the 
Naval Postgraduate School. An electronic method of obtaining con- 
tinuous graphical plots of probe current vs. applied probe potent- 
ial and probe current squared vs. applied probe potential was 
developed and successfully used. Preliminary experiments indicate 
this system may be extended to obtain the first and second deriva- 
tives of probe current with respect to applied probe potential 
allowing an accurate determination of space potential and" electron 
energy distribution functions. A gas diode filled with argon at 
5.9 microns of mercury was used in this investigation. Probe data 
was taken for twenty positions along the axial length of the tube 
using an anode-cathode spacing of 42 mm. Five different discharge 
currents were studied. A second run was made for ten positions 
and a discharge current of 50 ma at 6.9 microns of mercury to study 
distribution functions. Electron oscillations were detected which 
corresponded with the Langmuir plasma frequency predicted from the 
plasma density; the appearance of these oscillations corresponded 
with a decreasing electron density gradient. Peak intensity of 
oscillations occurred along a decreasing density gradient. 
Oscillations were also measured in the kilocycle range which are 
attributed to ion oscillations associated with beam electrons. 

Data is presented in the form of: raw data plots; space potential 

as a function of probe position; electron temperature as a function 
of probe position; electron density as a function of probe position; 
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electron density as a function of discharge current; and electron 
energy distribution functions. It is tentatively suggested that 
the unmodulated beam of charged particles, coming from the cathode, 
passing through and maintaining the plasma, suddenly becomes un- 
stable and gives rise to electron oscillations of the plasma in 
the presence of an associated decreasing plasma density gradient . 
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1. Introduction. 



The recent successes in harnessing the explosive fusion of 

the nuclei of the light elements and the even more recent attempts 

at space exploration has led to a renewed peak of interest in 

plasma physics in an effort to better understand the basic physical 

processes occurring in a highly-ionized, gaseous, conducting medium. 

Success in understanding these phenomena may ultimately lead 

to generation of fusion power through the medium of a plasma. 

Linhart (1) describes this possibility in the following way. 

It was appreciated that the energy output from fusion 
reactions depends critically on the kinetic energy of 
the colliding nuclei and that fusion outputs of 
practical interest depend upon one’s ability to produce 
temperatures of at least several million degrees Kelvin. 

If explosions are to be avoided, then the pressure of 
matter at this temperature must be balanced by external 
forces. This is within the power of our engineers only 
if the density of the nucLear fuel is less than the 
density of our atmosphere. The search for a mechanism 
of a controlled release of fusion energy became, there- 
fore, synonymous with the study of high temperature, 
low density plasmas. 

Since 99.9% of the matter in our universe is ionized and 
therefore composed of plasma, communication and electronic systems 
used in space exploration must be able to operate in, and propagate 
signals through, ionized regions of the outer atmosphere and space. 
Thus, an understanding of plasma processes will contribute strongly 
to the solution of these space problems. 
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Plasma physics is concerned with the behavior of systems of many 
free electrons and ionized atoms where the mutual Coulomb interactions 
cannot be ignored. In a restricted sense, such systems of particles 
consist of nearly equal numbers of positive ions and negative 
electrons Systems of this type are examples of a medium known as 
plasma which in many respects behaves differently from the solid, 
liquid and gaseous states of matter (1) . 

The average kinetic energy per particle determines the state in 
which matter exists. Mattel" can exist as a plasma, i.e., in its 
fourth state, when the kinetic energy per plasma particle exceeds 
the ionizing potential of atoms -- which is usually a few electron 
volts (1) . 

Plasma spans a broader energy band than any other state of 

matter; it occupies about 20 octaves of the kinetic energy spectrum . 

This wide kinetic energy band of the plasma state is the reason for 

so much common ground between plasma physics and many other fields of 

physics, such as the dynamics of single charged particles, or the 

physics of electrical discharges in gases (1) . Linhart says: 

iVlthough probably more than 99.9% of the matter in our universe 
is ionized and therefore in the plasma state, on our planet 
plasma has to be generated by special physical processes and 
special conditions. 

On earth the medium which most closely resembles an ideal plasma 
is a partially ionized gas. This fourth state of matter defined as a 
plasma is not very much at home on our dense, cold planet; therefore, 
special conditions must be realized in order to generate plasma-like 

conditions in the laboratory. These special conditions are in the 
realm of the physics of electrical discharges in gases (gaseous 
electronics) and hence, the close kinship between this discipline and 
plasma physics. 
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Plasma was originally defined by Langmuir and Tonks (2) 
in the following way: 

The word T plasma T will be used to designate that portion 
of an arc type discharge in which the densities of ions 
and electrons are high but substantially equal . It 
embraces the whole space not occupied by sheath. 

The positive columns of the glow and arc discharge are ex- 
amples of partially ionized plasmas. Plasma regions are highly 
conducting. In rare gas plasmas, the negative carriers are 
electrons. Ionization of the plasma is produced essentially by 
electron collisions with some photoionization also occurring. 

The positive ions, the electrons and the neutral gas atoms may or 
may not be in thermal equilibrium. The temperature of the positive 
ions is usually greater than the gas temperature, while the elect** 
ron temperature may be very high. Generally, positive ions and 
electrons may be assumed to have a Maxwellian distribution; 
however, variations from this distribution have been noticed 
(3, 4, 5, 6). 

Even though a plasma, with its ions and electrons moving 
at random according to their characteristic temperatures, is 
essentially neutral when considering a large volume, there must 
be high fields existing at points within the volume. Thus, at 
any point the field will vary with time depending upon the 
instantaneous random configurations of charged particles surround- 
ing the point . Fields of this type surrounding a point in a plasma 
may be called microfields and because of them, electrons may gain 
considerable amounts of energy in passing through a series of 
accelerating microfields. 
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Another property of an ionized gas which owes its origin to 
space charge phenomena and high effective conductivity is that of 
electrostatic plasma oscillations. These longitudinal osciLlations 
of the plasma are associated with the fact that a small momentary 
displacement of charges of one sign with respect to the other in the 
plasma will establish an electric field (microfield) in such a 
direction as to oppose the displacement. Inertial effects of the 
displaced charges then provide a mechanism for periodic osciLlations 
about the point of zero average charge displacement. The frequency 
of these oscillations is determined by the mass of the displaced 
charges and the f spring constant 1 of the electric field produced by 
the displacement (7) . 

A simple derivation of longitudinal plasma oscillations might 
be made as follows; 

We write Newton’s second law of motion, 

(1) d/dt(m dr/dt) = F =-eE. 

Now considering only wave propagation longitudinally, 
i.e., in the X direction, this reduces to 

(2) mX + eE x =0 

We also write Maxwell’s Curl equation 

(3) V x H = 4 TT J + . 

Since no external magnetic field is present, H = 0 

Now reducing the equation to its scalar component in the 
X direction, and denoting partial differentiation with respect 
to time by • , we have 

(4) 0 = - 4 IT neX + E x . 
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Then, integrating and setting the constant of integration 
equal to zero, since we have no external electric field present 



we obtain 

(5) E x = 4 TT neX 

Substituting (5) into (2) , we obtain 

(O') X + 4 t ne 2 X = 0 

m 

This equation gives for the frequency of oscillation 




Where: 

r is the vector displacement of the electron from 
the point of zero charge displacement 
e is the charge of the electron 
m is the mass of the electron 
F is the restoring force 
E is the electric field intensity 
H is magnetic field intensity 
J is the displacement current density 
t is time 

• is the partial derivative with respect to time . 

Plasma oscillations may be either electronic or ionic in 
nature, and the frequencies are determined by the respective 
electron and ion masses, hence, ion oscillations will be much 
lower in frequency than electron oscillations because of their much 
greater mass (7) . 

Plasma electron oscillations may transmit energy by a growing 
wave mechanism. A growing wave results from the excitation of 
plasma electron oscillations and the interaction of these oscillations 
with the primary beam. As the wave moves longitudinally through the 
plasma, it tends to accumulate the slowly drifting groups of oscill- 
ating electrons in advance of the wave and disperse the electron 

groups after passage of the wave. As this redistribution takes place, 
an exchange of energy between electrons and wave in a particular phase 
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relationship increases the energy in the wave. The wave amplitude 
grows exponentially reaching the maximum rate of growth when plasma 
frequency approaches the excitation frequency of the wave (8) . 

Energy of plasma electron oscillations may also be transmitted 
by fast-moving electrons traversing the oscillating region and having 
their velocities modulated periodically. This type of energy transfer 
obviously depends upon the relative velocities of the oscillating and 
non-oscillating groups of electrons (3). 

Experiments in plasma physics concerned with probe studies and 
plasma oscillations are usually carried out in gas-filled diodes or 
triodes utilizing discharges through rare gases or mercury vapor. 

A diode filled with Argon was used in this study. 

Typically, the voltage profile across such a tube from cathode 
to anode consists of a steep rise in potential near the cathode, 
followed by a slow gradual rise across the positive column of the 
discharge to the anode with a sharp, small rise or fall in the 
region near the anode. The electrons traversing the interelectrode 
distance gain their major potential rise in region adjacent to the 
cathode called the cathode Tr fall TT (9)^ The cathode is surrounded 
by a sheath of positive ions because of the greater mobility of the 
electrons compared with the ions. Electrons are accelerated as they 
move away from the cathode sheath producing more ionizations and 
attaining a Maxwellian distribution of velocities as they travel 
through the gas making collisions with gas atoms and each other. 

At the negative glow boundary, the entering electrons are composed 
of two or more groups -- fast electrons which have been produced at 
the cathode and have experienced few collisions in the cathode dark 
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space, and slow electrons which have made many Inelastic collisions 
in the dark space. The slow electrons have energies below the 
ionization potential of the gas, but at, or above, the excitation 
potential; therefore, they make exciting collisions with the gas 
molecules, causing light emission in the negative low region as the 
excited atoms emit quanta in returning to their ground states. 
Afterward, these slow electrons have such low energies that they 
then recombine with positive ions. The group of fast electrons are 
charge carriers which enter into the positive column part of the 
discharge. Some of these electrons may carry through to the anode 
depending upon the mean free path, but more likely, these electrons 
interact in the positive column, give up some of their energy and 
become randomized. The anode will then receive principally elect- 
rons coming from the positive column. 

This investigation was originally started at the United States 
Naval Postgraduate School in 1958 by Webb and Garner, continued by 
Flynn and Shoults in 1959, and by Culwell and Johnson in I960. The 
present work is principally a continuation of the work of Culwell 
and Johnson who redesigned and rebuilt the vacuum system and dis- 
charge tube, but were able to carry out only limited rudimentary 
measurements because of the press of time. 

The principal effort of this work was in (a) the design and 
development of a suitable electronic system with which to make 
intensive probe studies in the discharge tube and (b) then the 
making of such measurements. A unique system was designed yielding 
continuous graphical plots of conventional Langmuir data such as 
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probe current vs. applied probe potential and probe current 
squared vs. applied probe potential. This system shows promise 
of extension to yield accurate space potential determination and 
electron energy distribution functions. 

This report covers plasma oscillation and probe studies made 
in low pressure argon discharges. 
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2. Previous Work 



Investigation of gaseous discharge phenomena at low pressures 
was given its first real impetus in 1924 when Langmuir and Mott- 
Smith (10) reported a new method of making extensive studies of 
electrical discharges in gases at low pressures. Their technique 
consisted of placing a small auxiliary electrode (commonly referred 
to as a probe) in the discharge and studying the complete probe 
volt-ampere characteristic using the collector theory proposed 
by Langmuir (11) . 

Further studies by Langmuir (12) in 1925 disclosed that appar- 
ently there were three energy groups of electrons present in the 
discharge. He classified these groups as TT p rimary electrons TT 
(those electrons which keep nearly all the momentum gained in the 
cathode ion sheath) 5 TT secondary electrons TT (those electrons having 
a Maxwellian distribution of energies proportional to the energy of 
the primaries) 5 and a third group of TT ultimate electrons TT (those 
electrons having a Maxwellian distribution of energies independent 
of the primaries and much lower than the secondaries) . 

These observations indicated that there must be some type of 
energy transfer mechanism involved to permit certain electrons to 
gain energies in excess of the tube drop while others became energy 
deficient. Langmuir referred to this phenomenon as scattering. 

Floating probes had been used previously in analyzing discharges, 
however, they were of limited value because of the failure to appre- 
ciate sheath formation surrounding the probe and the consequent 
failure to measure space potential. 
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Dittmer (13) experimentally disproved the hypothesis that a 
new mechanism (such as the possible large increase in the Large t 
area of excited molecules) was the cause of scattering by showing 
that the scattering reached a constant value at some distance from 
the cathode. Dittmer attempted to find oscillations in an effort 
to validate the hypothesis that electrical oscillations in the dis- 
charge caused scattering. He was unsuccessful in all attempts at 
finding such oscillations. 

Tonks and Langmuir (2) developed a simple theory for electron 
and ion oscillations which was analogous to the electrons being im- 
bedded in a jelly-like substance composed of ions. The electrons 
have a natural frequency given by 



where n is the electron density per cubic centimeter, e the electronic 
charge, m e the mass of the electron. The ions obey a similar law with 
a natural frequency approaching 



where is the mass of the ion. 

Tentative verification of this theory was obtained by detecting 
oscillations in the 300 to 1000 megacycle range which corresponded 
reasonably well with the predicted frequencies. It was conjectured 
that additional oscillations in the 50 to 60 megacycle, 1.5 megacycle 
and lower ranges might correspond to oscillations of the primary 
beam electrons and ions, respectively. 
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Merrill and Webb (14) carried on extensive investigations in 
a mercury arc discharge at low pressures and discharge currents 
with a movable probe to determine a possible relation between the 
plasma oscillations and scattering as suggested by Tonks and 
Langmuir. They found that scattering occurred in several narrow 
regions and was associated with periodic oscillations of consid- 
erable magnitude which agreed with the Langmuir frequency within 
ten per cent. The fact that several regions were observed has 
been explained by Bohm and Gross to be a bunching effect of the 
electron beam. Merrill and Webb further noted that for each 
scattering region there was only one observable frequency and 
that the maximum intensity of the oscillations occurred at 
from one half to one millimeter from the scattering region on 
the anode side . 

Merrill and Webb concluded that the scattering was due to 

the oscillations receiving energy from the fast electrons , which 

then give some of the electrons excess energies, and that the 

observed difference in the positions of the oscillations and the 

* 

scattering was a consequence of the experimental technique . 



With the recent theory proposed by Allis and partially 
confirmed by Emeleus, one might find that an analysis of Merrill 
and Webb T s 600 complete volt-ampere characteristic curves would 
show that their observation between scattering and oscillation 
locations was of considerable physical significance. 
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Although no major breakthrough in understanding oscillations 
occurred in the next ten years, many important general observations 



were made. A few of these are stated below (15) . 

(a) The operating position on the current - voltage char- 
acteristic of the tube under investigation may affect 
the amplitude of oscillations observable in external 
circuits. It has been found that oscillations can 
usually be detected over wide sections of the 
characteristic for oxide-coated cathodes. 

(b) The position of the probe in the discharge markedly 
affects the amplitude of the observed oscillations; 
however, the observed frequency usually remains 
about the same. 

(c) The signal amplitude of the oscillations depends 
considerably on the electrode geometry and upon 
which electrodes are used for detection of 
oscillations . 

Emeleus (16) has pointed out the value of visual observations 
in electrical discharge studies using oxide-coated cathodes in ar- 
gon at a few microns of pressure with a low voltage discharge. A 
representative sketch of a gaseous discharge in which high frequen- 
cy oscillations are present may be shown as follows: 



-H / 




Fig. A. C, cathode; 1, cathode sheath; 2, unscattered beam 
section of discharge; 3, meniscus; 4,5, scattered beams; A, anode 
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The principal static fall in potential occurs in region 1; 
the rest of the discharge is essentially at anode potential. 
Between 1 and 2 represents the familiar Dittmer scattering region 
with dispersion occurring in or just after 3. It has further been 
noted that strong oscillations have only been detected to the 
right of 3 . 

Prior to 1950 little attention had been given to the effect 
of probe contamination on the validity of the volt-ampere probe 
characteristic. Easley (17) in probe studies of a mercury dis- 
charge buffered with argon showed, by drawing sufficient electron 
current to make the probe cherry-red (and thus outgas it by ohmic 
heating) before each reading that many of the non-linear logarith- 
mic plots previously obtained were due to a change in contact 
potential and thus did not indicate the presence of a non-Maxwell- 
ian distribution of electron velocities, but rather probe contam- 
ination . 

The results of Oleson and Found (18) have shown that non- 
linear logarithmic plots were obtained in discharges in which 
moving striations were present. The logarithmic plots obtained 
looked as though two Maxwellian distributions were present in 
the discharge. However, it had become fairly evident that in 
many cases (such as in discharges at pressures of the order of 
several microns of mercury) that a Maxwellian distribution of 
energies did not exist. Consequently, increased effort was 
directed toward new techniques of determining the actual energy 
distribution function rather than assuming it to be Maxwellian. 
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It is worthwhile to note that Langmuir considered the energy 
distribution function problem in his original theory on collectors, 
and though he did not employ the technique, he did show that ob- 
taining the second derivative of probe current with respect to 

* 

probe voltage would give the true energy distribution. 

Haigh and Smith (19) undertook a study of the energy distri- 
butions in the positive columns of low pressure (.005 mm Hg + 

.085 mm Hg) , hot cathode discharges in argon and air by a magnetic 
spectrum method. A continuation of this study by Haigh and 
McCulloch (20) obtained data on air, argon and nitrogen. Their 
results, using argon, are shown in Figs. B through E. These 
plots illustrate the types of distributions obtained in argon. 



Pressure 0.0825 mm Hg 
Tube Current - 15 ma 



Pressure 0.026 mm Hg 
Tube Current - 15 ma 




Energy (eV) 




Energy (eV) 
C. 



Fig. B. Fig. 

Strictly speaking, this is true for a cylindrical collector if 
the velocity distribution is isotropic. In cases in which one 
doubts the presence of an isotropic distribution, a spherical 
collector is required. This is shown in Langmuir’s original 
development . 
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It is interesting to note that for a given pressure, a decrease 
in current has the same effect as a decrease in pressure at constant 
current. Haigh and McCulloch hypothesized that the second peak in 
their distributions may be the result of the TT combined effects of the 
Ramsauer free paths and the tendency for a larger proportion of high 
velocity electrons as the pressure decreases. TT 

As far back as 1939 Druyvesteyn (21) had shown that the electron 
energy distribution function could be found from the expression 



/2m V 

N f (energy') =N0feV) = 2 / e P 

Ae"~ 1 

where N0(eV) is the number of particles per unit of energy per unit 
volume, and in this work, customarily put in units of reciprocal 
electron volts and reciprocal cubic centimeters, i.e., eV cm ; 
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A is the area of the probe, m is the mass of the electron, 

e is the electronic charge, ip is the probe current and Vp is 

the probe voltage. All arc expressed in consistent units. 

Despite this, no significant effort was expended to obtain 
energy distributions because of the inherent inaccuracies in- 
volved in the tedious graphical method of obtaining the second 
derivatives of probe current with respect to probe voltage. 

Sloane and MacGregor (22) developed an idea originally pro- 
posed by Emeleus for electrically obtaining the energy distribution. 
By putting a small alternating signal on the probe it was shown 

that the dc increase caused by the ac signal was directly propor- 

tional to the second derivative of probe current with respect to 
probe voltage. Though it appeared that this technique might yield 
the first real gain in obtaining distribution functions, Emeleus 
mentioned in Loeb (3) that "not too much confidence should be 

•k °k *-n 

placed in the distribution derived by this technique 
principally because of the many ways the plasma might react to the 
alternating potential. 

Recently, Boyd and Twiddy (6) used an improved version of the 
Sloane and Mac Gregor technique by amplitude modulating that ac 
signal thus permitting work in plasma having higher noise levels. 
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A typical energy distribution measurement by Boyd and Twiddy 
for a striated head in hydrogen shows at least two electron energy 
groups, as indicated in Fig. F. 




Fig. F. 

Boyd and Twiddy have also shown that for a non-isotropic, 
strongly directed discharge a cylindrical probe will give an energy 
distribution which will have negative values of energy between two 
major electron groups. In every case studied by Boyd and Twiddy, 
two well-separated groups of electrons (and sometimes a third group) 
were detected. Further, it has been suggested that the energy ex- 
change between the electron groups might be a much more rare phenom- 
enon than had been expected. 
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In 1956, Medicus (23) developed a simple technique for deter- 
mining the energy distribution function from the probe character- 
istic curves plotted automatically with an X - Y recorder. 

Medicus T method is a graphical one in which tangents are drawn to 
the probe curves and the differences between the tangents repre- 
sent the second derivatives. Fig. G below outlines the method 
pictorially # 




Fig . Go 



To obtain the best results the probe curve should be ampli- 
fied so that all regions of interest lie at about a 45° slope. 

The actual probe curves may be used directly only if the ion current 
is small compared to the electron current, or if the ion current 
varies linearly with voltage; otherwise, one must use the net 
electron current. 
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Medicus (24) obtained results which seem to further substanti- 
ate the suggestion of Boyd and Twiddy concerning the possible low 
energy exchange between slow and fast energy groups. Medicus showed 
that when a discharge consists of both a high energy , non-Maxwellian 
group of electrons and a low energy Maxwellian group that ,T the 
characteristic linear part of the corresponding probe curves, extra- 
polated to plasma potential, directly determines the random current 
density of the fast group.” Hie expression for this relationship 
is given by 4 i 



Fig. H shows how this is obtained graphically and how it is re- 
lated to the high energy electron distribution. 




e A 



i 



P 




Fig. H. 
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Applying the diffusion law to his result, Medicus found that 
the interaction between the slow and fast groups must be small and 
that cases may exist TT at low total plasma densities where the 
particle density of the slow group is always comparable to, or 
even much higher than, that of the fast group. 11 

Continuing the search to improve probe techniques considerable 
effort has been expended to develop efficient methods of eliminating 
probe contamination while obtaining the probe characteristic curves. 

Wehner and Medicus (25) , using an automatic X - Y plotter, 
were able to obtain a complete probe characteristic curve in sever- 
al seconds coming from a cherry-red probe condition. They claimed 
no deviation from a Maxwellian distribution over nearly 4 powers 
of 10 in probe current even in tubes with oxide-coated cathodes, 
providing the characteristic was recorded in several seconds after 
the probe had been thoroughly cleaned as described above. 

White (26) , while making measurements in a hollow-cathode 
glow discharge, found that excessive material was sputtering from 
the cathode onto the probe. It was found, that by keeping the 
probe at 30 volts below floating potential and by pulsing it 
every 0.1 seconds for 0.004 seconds, it was possible to maintain 
a clean probe by ion bombardment and hence, obtain a reliable 
p^obe characteristic . 

‘k 

' This would imply that the particle densities of the different 
energy groups would have discrete values at any point in the 
discharge tube for a given set of parameters. 
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Several months later Waymouth (27) reported a similar 



pulse technique for obtaining "uncontaminated" probe 
characteristics . 

Some results of electron current vs. probe voltage For 
discharges at pressures of six microns of mercury plus one 
mm of neon are shown in Fig. I, below, taken froir Waymouth’ s 




Curve 

A 

B 

C 

D 

E 

F 



Distance 
From Cathode 
0.25 cm 
0.75 cm 
1.50 cm 

2.00 cm 

3.00 cm 
4.25 cm 



An attempt to obtain energy distribution functions by the 
method of Medicus gave poor results because of the small size of 
the oscilloscope photographs. Standard scatter of data is shown 
for a probe 0.25 cm from the cathode. 
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‘A' 

It was found that this pulse technique represents at 
least a 600-fold reduction in the net rate of probe contam- 
ination in gas discharges. It was also shown that the energy 
distribution of electrons in the negative glow consists of 
two groups, a Maxwellian group and a high energy group with a 
potential approximately equal to the cathode fall. 



It is important that in the zeal to overcome one problem 
it is not replaced with a more complex one. For these measure- 
ments to have meaning the ion sheath must be in equilibrium; 
therefore 5 it is clear that there is some upper limit on the 
rise time of the pulse so that the sheath thickness is essenti- 
ally in constant equilibrium, i„e., where the distance the ions 
can travel during a pulse equals the product of the pulse dura- 
tion, the ion mobility and the sheath field. 
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Waymouth (28) has pointed out an important new 

diagnostic tool in making a rapid evaluation of the 

principal electron energy groups directly from thp linear 

i vs V„ probe characteristic curves. His technique is 
p t J 1 

essentially this . 

N0(eV) ~ d^i p /dVp , 

therefore, whenever the linear i p vs V p curve has a constant, 
slope 

N0(eV) = 0, 

i.e., there are no electrons of corresponding energy in 
the energy distribution spectrum. 

An example is shown in the following plots; 




Fig. K 
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A modern radio frequency technique of determining the 

plasma resonant frequency in a gaseous discharge has been 

developed by Takayama, Ihegumi, and Miyozaki (29) » They 

discovered that if an rf signal is superposed on a probe, 

which is maintained negative at all times relative to space 

potential, then a resonant current increase appears in the 

circuit at the Langmuir frequency, f = I ne^ 

y1Tm e 

The relation between the resonance and the frequency of 
imposed rf voltage is as shown in Fig. L, below, for a 
homogeneous mercury plasma at 0.1 mm Hg pressure. 





Takayama, et al., suggest that this is a beam-plasma 

type of interaction and that the resonance effect will not 

be discernible when the energy supplied by the rf field is 

comparable to the thermal energy of the electrons. They 

point out that an increase in pressure will likewise cause a 

decrease in resonance since the height of the resonance peak 

also depends upon the collision frequency of the electrons with 
atoms. 
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In 1957, Allis (30) suggested that as an electron beam 
travels through a region of increasing plasma density, electron 
oscillations are being generated in a region in which they can- 
not propagate in the absence of the beam. On the other hand, 
an electron beam travelling through a plasma of decreasing 
density might generate oscillations which could propagate and 
outrun the beam. 

Allis further pointed out that if a monoenerge tic beam of 
electrons is going through a plasma then the average beam velocity 
is reduced by the growing amplitude of the oscillations, the 
energy being transferred into field energy in the wave and in- 
creased kinetic energy of the random plasma electron group. 

Allis states, "growth occurs mainly below the plasma frequency" 
so that the free plasma cannot transmit these waves, (see Fig. M) . 
He further mentions that if the beam should peter out through 
scattering or slowing down, the plasma waves must decay expon- 
entially along with the beam. 



Amplification 

Factor 




Fig. M. Amplification factor for an electron 
beam in a plasma (30) . 
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A recent investigation of electron beam interaction with a 
plasma was conducted by a group of Russian workers (31) . In their 
paper they point out that theoretical investigations have shown 
that the state of an unmodulated beam of electrons moving through 
a homogeneous plasma is unstable. Density and velocity fluctuations 
arising in the beam or the plasma, may excite longtudinal electro- 
magnetic waves (plasma oscillations) of growing amplitude in the 
plasma. They note that the initial fluctuations may be caused by 
external forces. 

These workers state that the observed maximum excited or 
amplified frequency is determined by the following relation: 



^o = 




’ for V o 3> s 



Where , 

^o is the observed frequency in radians 
cOp is the Lnagmuir plasma frequency 
is the beam velocity 

S is the average thermal plasma velocity 
They also point out that a large velocity spread in the beam 
may remove the instability. If the regularly moving beam consists 
of ions, the excitation is less, not by a factor m e /m-£, but only 
(m e /m^) . Therefore, the interaction mechanism considered is 

also important if an ion beam moves through a plasma. 
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Plasma Instabilities may be either convective or absolute. 
Convective instabilities tend to move with the flow and the 
oscillations they produce amplify or diminish with time. Ab- 
solute instabilities produce oscillations which grow contin- 
uously with time. 

Emeleus and Garscadden (32) have very recently done some 
work in exploring plasma oscillations with a transverse electron 
beam observing the resultant energy spread of the transverse beam 
when it is in, or has traversed, a section of the main discharge 
where scattering or oscillations can be detected. 

A typical result, for a discharge of 1.5 microns pressure 
in mercury vapor with a discharge current of 27 ma (8 mm diameter 
cathode) with an electrode separation of 42 mm and a meniscus 
forming 10 mm from the cathode, is shown below. 




Fig. N. Energy spread for a transverse electron 
beam in a meniscus-type discharge. 
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The conclusion is made that "with the discharge described 

* * * The results appear to be compatible with considerable 

randomization and possible break-up of the main beam by the 

oscillation" as the beam interacts with the plasma in the 
* 

meniscus region. 

Very recent investigations by Mahoffey, McCul.Lagh, Garscadden 
and Emeleus (33) are pertinent in correlating growth or decay of 
oscillations with plasma density gradients. 

Their work was also done in a mercury discharge at approxi- 
mately one micron of pressure using an eight mm oxide-coated 
cathode with a tube current of 27 ma and a tube voltage of about 
20 volts. Two main sets of measurements were made with emphasis 
on (a) frequencies 3 intensities, and energy spread of the electron 

beam (found from the break in the i vs. V n curve at which fast 

p r 

electrons began to be received) . (b) frequency measurements, 

electron temperatures, and electron concentrations (found from 
ip vs Vp curves near space potential) . 

They conclude that: 

Some correlation has been found between growth or 
decay of oscillations and motion of the electron 
beam down or up a plasma concentration gradient, 
which appears to be in general accord with the 
plasma dispersion theory. 



*cf . our results, Section 5 of this report, for a beam generated 
in the discharge tube interacting with the plasma. 
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3 . Equipment . 



The equipment used can be divided into four main component 
groupings as follows: 

(1) Gaseous discharge tube in which a glow discharge could 
be continuously maintained. 

(2) Vacuum system for evacuation and filling of the tube. 

(3) Frequency detection and measuring equipment. 

(4) dc electrical circuit for supplying power and for 
controlling and measuring discharge characteristics. 

Figs. 1 through 3 are photographs of the general experi- 
mental setup 

a. Discharge Tube . Figs. 4a and 4b are photographs and Fig. 

5 is a detailed drawing of the tube (diode) . This is the same 
tube used by Culwell and Johnson in their work during 1960. No 
changes were made to the tube. The electrodes and probes were 
attached to soft iron cylindrical cores sealed in glass which 
made possible the variation of interelectrode distances and 
probe positions by moving the iron cores with externally 
applied magnetic fields. 

The probe used in this work was a cylindrical probe made 
of .008 inch diameter tungsten wire enclosed in a glass sheath 
up to 3.8 millimeters from the tip. 

The ring probe shown in Fig. 5 was not used during the 
experiment. It was damaged during electrode movement, breaking 
loose from its mounting and falling to the bottom of the discharge 
tube . 
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The cathode was a Philips impregnated TT L TT cathode with a 
two centimeter diameter face indirectly heated by a tungsten 
element drawing seven amperes dc . 

The anode was a three centimeter diameter nickel disc. 

This tube was originally designed with high purity of the 
gas being a primary design criterion. As such, the coiled wire 
leads from the electrodes and the probe to the external connect- 
ions were unavoidable. The inductive effects of these leads may 
be considerable at electron oscillation frequencies and tend to 
mitigate against making large amplitude high frequency measure- 
ments through them. 

b. Vacuum System . The vacuum system built by Culwell and 
Johnson (8) was not changed. It remained operative during the 
entire course of this work without leaking. Pressures of the 
order of 10” mm of mercury were readily attained. All valves 
in the system were glass stopcocks lubricated with Type N Apiezon 
grease. The diffusion pump was a three-stage air-cooled pump 
using Octoil-S. The major part of the system was baked out with 
a portable oven (see Fig. 1 ) . Heating tapes were used to bake 
out external glass tubing. An induction heater was used to out- 
gas the anode. Outgassing of the Langmuir probe was accomplished 
by ohmic heating while the tube was operating. 



* 



Manufactured by Semicon of California, 



Inc . 
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Low pressures were measured with a Veeco Type RG 3 L A 
ionization gage. Pressures in the operating range of one to 
ten microns were measured with a Langmuir viscosity gage (35) „ 

Tube voltage readings were also used as an indication of tube 
pressures after the voltage-pressure relation had been determined. 
Tube voltage is a function only of tube pressure if other para- 
meters such as anode-cathode spacing and discharge current remain 
constant; therefore, a given pressure could be readily duplicated 
by evacuating or filling the tube until the correct voltage was 
obtained. Thus, when fresh argon was desired in the tube, flush- 
ing and refilling was quickly and accurately accomplished without 
relying on the slower pressure measurement process using the 
viscosity gage. 

c. Frequency Detection Equipment . A cavity wavemeter. Test 
Set, Model AN/UPM - 2, with a frequency range of 80 to 1220 Mcs 
was used to detect electron oscillation frequencies by coupling 
with the wire probe at anode potential. The dc signal from the 
cavity wavemeter detector crystal was fed into a Hewlett-Packard 

Model 425 A dc micro volt-ammeter to obtain detectable signals. 

-12 

Currents as small as 10 amperes could be detected with this 
instrument . 

A Hewlett-Packard UHF signal generator, Model 612 A, with 
a frequency range of 450 to 1230 Mcs was also tried utilizing the 
technique of parametric amplification suggested by K. Takayama, 
et. al. (29) which essentially consisted of putting the rf 
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potential on the probe held at space potential while the frequency 
of the generator was varied in an attempt to cause resonance at 
the electron oscillation frequency. 

Ion oscillations were detected with a Tektronix Type 581 
oscilloscope and also simultaneously counted with a Hewlett-Packard 
Electronic Counter, Model 521 A. The ion signals were amplified 
by a Hewlett-Packard 40 db amplifier. Model 450 A, to feed count- 
able signals to the counter. 

d. dc Circuitry Components . Fig. 6 is a schematic block dia- 
gram showing the general dc circuitry employed. 

Filament voltage was supplied by storage batteries regulated 
by a charging circuit during operation. Seven amperes dc was 
sufficient to heat the Philips L type cathode to 1050° C. Cathode 
temperature was measured on the face with a Leeds and Northrup 
optical pyrometer. Cathode temperature was maintained within + 

5° C during operation. 

Electrode spacing was measured with a locally manufactured 
cathe tome ter which allowed electrode positioning within + 0.25 mm. 

The discharge tube was powered by a Kepco Labs Regulated 
High Voltage Power Supply, Model 1250 B, having a voltage range 
from 0 to 1000 volts and a current output from 0 to 500 milliam- 
peres . 

An F. W. Moseley Autograf X - Y Recorder, Model 2, was used 
in the plate circuit to obtain continuous plots of the volt - 
ampere tube characteristic simultaneously with the taking of probe 
data, thus the effect on discharge current of drawing probe 
current could be monitored and recorded. 



32 



lateral attempts at taking probe data were designed around 
oscillographic recording of the probe characteristic curves. 

Spurious high frequency signals in the range 0 S to 5 Mcs of un- 
known origin interfered with these attempts despite coaxial 
shielding, etc. The oscillographic method was therefore abandoned 
as a primary means of measurement and a system designed to mini- 

c 

mize ac interference was constructed. A description of this system 
follows . 

The probe voltage was supplied by a Hewlett - Packard Low 
Frequency Function Generator, Model 202 A, set for triangle wave 
generation to obtain a voltage source varying linearly with time. 
This signal was then amplified and biased to obtain the desired 
positive and negative voltage excursions relative to the anode, 

A Donner Analog Computer, Model 30, with Problem Board, Model 
30-6, was used to amplify the voltage drop across the variable 
resistor in the probe circuit, thus, producing a voltage proportion- 
al to probe current in order to ensure a voltage range suitable for 
best operation of thq function multiplier and X - Y recorder. Ten 
operational amplifiers were available with this computer. The 
analog computer was also used as a source of five separate regulat- 
ed dc power supplies from its initial conditions panel. The analog 
computer amplifiers operate in the range -100 to +100 volts with 
long time stability better than 0.5% full scale value. Initial 
condition voltage sources operate from 0 to 100 volts. Computer 
amplifiers have a maximum drift characteristic of 4 millivolts 
per hour, with a short term random drift less than + 2 mv. 
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Balancing of amplifiers is easily accomplished and holds drift 
error to negligible levels over short time intervals for high 
signal to drift ratios. Potentiometers on the problem board 
allowed precise gain control on amplified signals. 

A Donner Electronic Function Multiplier, Model 33, was used 
to square the values of voltage proportional to probe current. 

This equipment would accept inputs of + 100 volts in either or both 
input channels (X, Y) and produced an output voltage of -0.01 XY. 

An F. W. Moseley Autograf X - Y Recorder, Model 2S, was used 
to make continuous plots of probe data. Probe voltage with re- 
spect to the anode was applied to the horizontal axis which has 
a range from 0-7 .5 millivolts to 0-150 volts full scale. The 
voltages proportional to probe current, or probe current squared, 
were applied to the vertical axis which has a range 0-5 millivolts 
to 0 - 100 volts full scale. Input impedances of both axes are of 
the order of 2 megohms. Sweep rate was controlled by the low 
frequency function generator. It was generally kept from 0.05 
cps to 0.2 cps to avoid inertial error in the X - Y recorder. 

This error is negligible at these sweep rates. The Autograf is 
accurate to 0.25% of full scale value on all ranges when calibrated. 
Calibration drift is negligible. 

Both the output of the function multiplier and input of the 
X - Y recorder have low pass filters keeping ac interference to 
a minimum. 
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A Lambda Electronic Corporation Regulated Power Supply, 
Model 71, was used as a probe voltage supply when a non time- 
varying source was desired. 

A Donner Variable Base Function Generator, Model 3750, was 
also tried in an attempt to automatically make plots of the 
natural logarithm of probe current vs. applied probe potential. 
Preliminary attempts were unsuccessful because of difficulty in 
programming a suitable logarithmic function. This technique was 
temporarily abandoned because of the press of time. It does 
appear quite feasible, nevertheless, and pursuit along this line 

k 

might prove to be quite rewarding. 



* 

F. W. Moseley Co., manufactures a logarithmic attachment for 
the X - Y recorder which enables semi-logarithmic plots to be 
made without need of a function generator. They also manufacture 
a curve follower attachment to the X - Y recorder which allows 
functions to be drawn on paper with magnetic ink; the pick-up 
coil will then feed out the function desired as the independent 
voltage is driven onto the horizontal axis. A system such as ours, 
utilizing this technique, would require two recorders -- one to 
record, one to generate functions both being driven in time 
synchronism. 
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4. Procedure . 



Prior to making a data run, the vacuum system was evacuated 
to as low a level as possible — in the range of 10~ 4 mm of mercury. 
Then the system was baked out with the oven at 450° C for a period 
of between twenty four and thirty six hours. During this time the 
tubing and liquid air traps external to the oven was wrapped in 
heating tape and heated to 350° C. Following bake out, the cathode 
was heated to operating temperature with the filament and the anode 

•k 

was outgassed with the induction heater. The system was then 

8 -9 

allowed to pump down to between 5 x 10 and 10 mm of mercury. 

This normally took 24 hours after completion of bake out. The 
system was then filled with approximately 2 mm of spectroscopic- 
ally pure argon. This filling procedure was repeated several 

times in order to flush the tube with the inert gas. With argon 

-8 

in the system, the tube could be evacuated down to 5 x 10 mm 
of mercury in about one and one half hours. The final pressure 
was obtained by discharging the tube, measuring the pressure with 
the Langmuir viscosity gage, and correlating it with the plate 
voltage; then, evacuating the tube in a small increment of pressure 9 
while the tube was still discharging, the process of measuring 
pressure and correlating with tube voltage was repeated until 
pressure - tube voltage calibration had been achieved. The final 
pressure was then readily obtained by evacuating until the proper 
plate voltage was obtained. 

k 

This bake out and purification procedure is tedious, but necess- 
ary. For example, failure to be careful about extreme cleanliness 
results in a pulsating discharge with a frequency of one to 25 cps 
and an extremely unstable tube characteristic (i.e., the tube ’’runs 
away 1 ’ at very low discharge currents) . 
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An arbitrary interelec trode separation distance of about 42 
mm was selected so that probe measurements could be made with the 
probe at increments of two mm along the axial length of the tube 
between electrodes. Thus, about twenty different positions between 
anode and cathode were studied at approximately equal intervals. 

Electrode positions were calculated and set at the catheto- 
meter as follows: The cathetometer was set at its reference value 

for the cathode. The cathode was then moved using the magnetic 
coil over the glass surrounding the soft iron magnet within the 
tube until the cathode was lined up with the reference mark on the 
cathetometer. The cathetometer was then traversed to the probe 
position where its value was recorded (since the probe position 
remained fixed, this position was used to estimate the statistical 
accuracy of the cathetometer measuring system) . The cathetometer 
was then traversed to the anode position and the anode moved to 
coincide with the reference mark as described for the cathode . 

All cathetometer readings were made going in the same direction 
to minimize backlash errors. 

A tube volt-ampere characteristic trace was then made and 
recorded on the plate circuit X - Y recorder to determine the 
operating range of tube currents. When the gas in the tube was 
very clean the tube remained stable up to 200 milliamperes dis- 
charge current. Over this range, then, a group of representative 
discharge currents were selected to study which gave suitable 
values of probe current (i.e., probe currents which did not pro- 
duce voltages in excess of 100 volts to the operational amplifiers 
of the analog computer') . The discharge currents selected were 
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25, 50, 75, 100 and 125 milliamperes . All of these currents were 
attainable using the voltage control on the Kepco Power Supply, 
the plate circuit resistance remaining fixed. 

Proceeding next to the probe circuit, calibration procedures 
were employed to minimize all electronic errors. The amplifiers 
of the analog computer were balanced to minimize drift. The 
function multiplier was balanced for drift and accurate suppression 
of cross products. A potentiometer on the problem board was used 
to calibrate the maximum peak-to-peak excursion of the triangle 
probe voltage supply to an absolute magnitude of 150 volts as 
traced on the horizontal axis of the X - Y recorder. This voltage 
was then biased with an initial condition voltage source (set with 
a Donner Null VTVM) to give the maximum positive and negative ex- 
cursions desired. Those used were +20 volts to -130 volts with re- 
spect to the anode. 

Calibration of the vertical axis voltage corresponding to the 
probe current was accomplished with a problem board potentiometer 
after the signal had been amplified in the two nominal gain- ten 
stages. This potentiometer was then adjusted to give the same 
magnitude at penpoint on the X - Y recorder as the ammeter reading 
probe current. This had to be accomplished dynamically, but could 
be done to quite good precision (as well as one could read the 
meter) at very slow sweep speed. This was commonly done at the 
maximum ammeter and pen excursion when both momentarily came to 
rest. Inertial effects were considered to be negligible. 
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Once these calibration procedures had been made, data runs 
were started. For each position of the probe a linear probe 
characteristic curve was made for each of the five tube discharge 
currents cited previously. Blowup traces of these curves were 
also made of the ion (negative) portions of the trace. The paper 
was then changed in the X - Y recorder and the current squared vs 
applied probe potential plots were made for each of the discharge 
currents. Different vertical scales were used for electron and 
ion regions of the curves corresponding to those used for the 
linear traces. Simultaneous with these measurements, the plate 
circuit volt-ampere characteristic was being monitored and traced 
giving a record of its changes during the probe characteristic run. 

Prior to each new position of the probe the calibration 
procedures cited above were employed, except that analog computer 
amplifiers were only balanced approximately hourly. 

A run normally took twenty to twenty five minutes to make -- 
about five minutes to trace the curves, and about fifteen to twenty 
minutes to reposition electrodes and calibrate equipment. Twenty 
such positions between anode and cathode were studied using this 
procedure during a twelve hour period yielding four hundred 
separate probe characteristic traces. 

Reproducibility of curves is quite good within the limit of 
electronic errors which will be commented upon later in this report. 
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A subsequent, more limited, data run was made for ten equi- 
spaced probe positions (electrode separation 40 mm, argon pressure 
6.9 microns of mercury, discharge current 50 ma) , for the purpose 
of obtaining second derivatives of probe current with respect to 
probe voltage. Two simple RC differentiating circuits in cascade 
were used to do this on the probe current signal, (see Fig. 6) . 

Thse derivative curves were used then to obtain electron 
energy distribution functions. 

On this data run, ion oscillations were measured by 
oscilloscopic photographs and electronic counter measurements. 
Electron oscillations were measured with the cavity wavemeter 
coupled to the dc micro-ammeter . 

During this run great care was taken to evaluate the effect 
of probe contamination by retracing a curve after various time 
intervals and noting trace deviations from the original; then 
after cleaning the probe by ohmic heating, the curve was again 
retraced to note the tendency to return toward the original 
trace. The results of these tests were quite conclusive and 
are stated later in Section 5, Results. 
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5. Results. 



Plots of two representative probe positions from the first 

run (consisting of twenty probe positions and five discharge 

currents) are shown in Figs. 7 through 13 for the following 

functional relationships: 

2 

ip vs V p , ip vs In i p vs V p , 1^ vs V^, where the sub- 

script p denotes the probe parameter, subscript + denotes the ion 
region of the probe characteristic curve, subscript - denotes the 
electron region of the probe characteristic curve, and subscript 
d denotes the tube discharge (or plate) parameter. 

A summary of three important plasma parameters for all 
twenty probe positions between anode and cathode are shown in 
plots as indicated below: 

(a) n (electron density) vs X (cathode - probe distance) 

2 

obtained by using the slope of i p vs V p curves (Langmuir Theory) 



V 



is shown in Fig. 14. 

(b) V gp (space potential) vs X, obtained from In ip vs , p 
plots at the first positive voltage deviation from linearity is 
shown in Fig. 15. 

(c) n._ vs J (discharge current density) for ten probe 
positions spaced approximately 4 mm apart to correspond with the 
positions used during the second data run is shown in Fig. 16 . 

These data are perhaps capable of being correlated with the re- 
sults of the second data run, which was made for only one discharge 
current (50 ma) , despite the fact that the two runs were made at 
slightly different pressures. 
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A second data run was made at 50 ma discharge current for 
10 probe positions spaced at 4 mm intervals with an anode - cathode 
separation distance of 40 mm. The gas pressure was 6.9 microns of 
Hg. Results of this run are indicated below: 

(a) Plots of ip vs Vp, dip/dVp , d 2 i p /dV 2 are shown in 
Figs. 17 a, b through 18 a, b for two representative probe posi- 
tions . 

(b) Electron energy distribution functions for all ten probe 
positions are shown in Figs. 19 through 28. These functions were 
obtained from the second derivative of probe current with respect 
to probe voltage according to the relation 



N0(eV) 




Where N0(eV) is the number of electrons per electron volt of energy 
per cubic centimeter, A is the area of the probe, m is the mass of 
the electron, e is the electronic charge. For the geometry of the 
probe used in this work, the function reduces to 



N0(eV) 



13 % 

= 1.722 x 10 V 

P 





-3 

cm 



Two things are to be noted concerning these derivatives and 
associated energy distribution functions. First, the differentiat- 
ing circuit was not perfect — in particular it appears that the 
capacitor was not unloading quite fast enough to give true values 
when the derivative was decreasing. This error was relatively 
small and was arbitrarily adjusted to make it fit the first 
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derivative. Second, the X - Y plotter developed some mechanical 
troubles during this run and gave a magnitude for the second 
derivative of approximately 1/3 of the true value. Approximate 
calibration of the derivatives was achieved by using another 
machine which was reading true and also numerically differentiating 
(well realizing the errors inherent in this process) the i vs V 

Jr Jr 

curve to establish the scale factor. The second derivatives and 

the distribution functions shown in Figs . 19 through 28 are shown 

in their original scale as they came from the machine. The 

correction factor is applied to the constant term, i.e. K corrGC ^ on 
13 ^ ^2^ 

x 1.722 x 10 V p p This correction factor may still be 

dV^ 



in error by perhaps ten percent; however, these results are shown 
because the technique appears valid and the information the curves 
yield appears to be a valuable diagnostic tool. The actual magni- 
tudes really are not critical for diagnostic purposes; shapes and 
changes in shape in going from one probe location to another are 
important because these factors show how the energy of the particles 
is being transferred. 

(c) A plot of n_ vs X is shown in Fig. 29. These values 
were obtained by integrating the energy distribution functions 0 By 
this method the density could be divided between the high and low 
energy groups and changes in density groupings noted as the probe 
position was changed. The integration was carried out numerically 
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Results were checked with a plan!- 



& 

using the trapezoidal rule, 
meter and agree within three per cent. Fig . 47 shows a typical 

distribution function, and the technique employed to integrate it 
and obtain the various densities. 

(d) A plot of n_ vs X is shown in Fig. 30 <, These data were 
obtained by the conventional Langmuir method of using the slope of 
ip vs Vp plots. Comparison of Fig. 29 and Fig. 30 shows that both 
plots have the same general shapes, and essentially the same peaks 
and gradients are obtained by both methods. The absolute magni- 
tudes vary by as much as 30%. This may be attributable to the 
scale factor assumed for the second derivative method; but, possibly 
part of the difference (perhaps ten per cent) may exist in comput- 
ing the slopes by the Langmuir method since correction terms were 
not used, and only a relatively few points were used to determine 
these slopes. Since the main purpose in obtaining these values 
was to check the approximate magnitude of those values obtained from 
the distribution functions, great significance is not attached to 
these small differences. Even though absolute magnitudes of the 
densities computed from the distribution functions are not certain, 
the magnitudes relative to each other are thought to be quite good 
and this fact is used later to correlate density with the predicted 
Langmuir frequency and to make a statistical estimate of the 

ft 

The trapezoidal rule gives an approximate value to the area under 
an irregular curve, i.e.. Area = ^( y /2 + y 2 + . . . + y n ± + y n /2) . 

In our integration we used a 2^ of two electron volts giving a fair 
approximation to the integral. Using smaller increments gives an 
increasingly better result. Simpson T s rule is more accurate 0 We used 
the trapezoidal rule because it is suitable for making a running tabu 
lar integral. Simpson’s rule is not conducive to this treatment. 



standard deviation for the values of density obtained, 

(e) A plot of relative electron density vs X is shown in 
Fig. 31. These data were obtained from In i vs V curves by 
using the current at space potential and the electron temperature. 

(f) A plot of relative electron density of high and inter- 
mediate energy beam electrons is shown in Fig. 32. These data 
were obtained by Waymouth T s technique of analysis of i^ vs Vp 
probe characteristic curves to obtain the energy of the group* 
Medicus T technique, based upon a diffusion law which is dependent 
upon pressure, was used for relative magnitudes . Our pressure 
range was significantly different (lower) from that of Medieus ? 
development; hence, the results are good only to an order of 
magnitude . 

mm 

(g) A plot of eV (average electron energy) vs X is shown in 
Fig. 33. These data were obtained from the average energy of the 
two electron groups. The average energy of each group being 
obtained from the relation 

eV 0 



eV = 



/ 

eV. 



v 3/2 d 2 i 



d V„ 



d(eV) 



eV, 



eV P 



d 2 i 



■*? 



d(eV) 



These values are shown on the distribution functions in Figs. 19 
through 28. 
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(h) A plot of (electron temperature) vs X is shown 
in Fig. 35. These data were obtained from the slope of the 
linear portions of the In i vs V curves after the current 

lr ir 

attributable to the high energy electrons had been subtracted out. 

(i) A plot of T e vs X for both high and low energy electron 
groups is shown in Fig. 36. These data were obtained from the 
average energy of the two electron groups by using the relation 
eV = 3/2 kT e . This is strictly true only if the distribution is 
Maxwellian. The distributions obtained were close to this for the 
low energy groups and the relation was used for the high energy 
groups as well for want of a better relation to use. This plot is 
essentially the same as Fig. 33, since both were obtained from 
the same data. The ordinates differ only by a constant. 

(j) A plot of V S p vs X is shown in Fig. 37. These data 
were obtained by taking the maximum of the first derivative for 
positive values of V as being space potential. This corresponds 
to the inflection point of the i^ vs Vp curve, or the first devia- 
tion from linearity on the In ip vs Vp plot. 

(k) A plot of V vs X is shown in Fig. 38. These data were 
obtained by taking the first deviation from linearity on the In ip 
vs Vp plot as being space potential. Agreement between this figure 
and Fig. 37 is quite good. More faith is placed in the values ob- 
tained from the first derivative because these values are easier to 
read precisely than the deviation from linearity on the logarithmic 
plot . 

(l) A plot of E (electric field strength) vs X is shown in 
Fig. 39. These data were obtained by taking the changes in spdce 

potential as a function of distance (X) from Fig. 37 
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Plasma oscillations in two frequency ranges were found. 
Kilomegacycle signals detected are believed to be electron 
oscillations and kilocycle signals are believed to be ion 
oscillations . 

High frequency oscillations of 472 Mcs and 520 Mcs were 
found at positions 20 mm and 24 mm, respectively, from the cathode 
for the 6.9 micron, 50 ma discharge. Computing the density by use 
of Langmuir's theory relating the density to the slope of i^ vs 
Vp curves, we obtain the plot of total density shown in Fig. 30a, 
below. The densities corresponding to the observed frequency, 

i, 

calculated by using Langmuir's relation f = 8980 n 2 , are shown 
as dashed lines on this plot. 



m 

E 




n l- 

n 2_ 



Fig. 30a. Electron density as a function of cathode - probe 

distance. Data obtained from slopes of \_ c vs \L. plots. 

P p ' 
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Neither of these densities, corresponding to the observed 
frequencies, correspond with measured points on the totaJ density 
plot, but instead, represent plasma oscillation frequencies lower 
than would be obtained from total density values. 

However, if we make use of the density plot obtained from the 
energy distribution functions and also use Waymouth T s technique to 
obtain the distribution of high energy beam electrons, we obtain 
the plot of densities shown in Fig. 29a, below. By using these 
techniques we are able to break our electron concentrations down 
into groups which are: high energy beam electrons (obtained by 

Waymouth’s technique) (see p. 23, Section 2, Previous Work), 
intermediate energy beam electrons (obtained from energy distri- 
bution functions) and random plasma electrons (obtained from 
energy distribution functions) . 
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Fig. 29 a. • Total electron density (from distribution functions) 
O Plasma Density (from distribution functions) . 

, X Intermediate energy beam electrons (from 
distribution functions) 

— High energy beam electrons (Waymouth T s technique) . 
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If we apply the relationship cited by Kharchenko, et aJ . 
(31) (described on page 26 , Section 2, Previous Work) . 

CD *o = f p 




where S is the thermal plasma velocity and v q is the beam 
electron velocity and put this in terms of energy, we have 




where eVp is the average thermal plasma energy and eV^ is the 
average energy of the beam electrons. 

We know f Q , eV^, and eV^ experimentally. Let us now calculate 
the equivalent density corresponding to fp and see where this fits 
on the density plot. Fig. 29a. 

Consider first f G = 520 Mcs . 



(i) 




5.5/22 



^^intermediate 



= 0.25 



(4) and, 



f 

P 



% 

(1 - 0.25) 520 Mcs 

451 Mcs. 



: . v 

Realizing this was developed for v S 
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The equivalent density is 



( 5 ) 




451 x 10 
~ 8U5U 



6 



2 




Similarly, for f Q =472 Mcs, 



(6) eV, 



= 5/50 = 0.1 



eV, 

b high 

% 

(7) and, f = (1 - 0.1) 472 Mcs 

Ir 



f 



447 Mcs. 



P 



The equivalent density is 



(8) n 2_ - 2.47 x 10^ particles/cm^ 



Referring to Fig. 29a, densities n^ and ri 2 (and their 
corresponding frequencies) are indicated by dashed lines. It is 
apparent that these densities correspond with the curve repre- 
senting the plasma density, not total density . 

In obtaining the data for Fig. 29a,we were beset by a 
dilemma. We have previously mentioned that the densities ob- 
tained by the second derivative method were in error by some 
constant factor caused by improper machine calibration. However, 
we are confident that the relative values of these densities arc 
rather good with respect to each other. Now the absolute values 
of densities obtained by the Langmuir method are probably accurate 
to about 10%. In order to use our density data obtained by the 
second derivative method to best advantage we have determined a 
scale factor to be applied to these densities (obtained by the 
second derivative method) by taking ratios of each corresponding 
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ordinate of Fig. 29 (second derivative data) and Fig. 30 
(Langmuir data) . We have used the average of these ratios 
to correct both curves to the same absolute magnitude. 

Then, by statistically analyzing the ten density ratios so 
obtained, we are able to compute a number which gives a certain 
measure of confidence in the values of density plotted in Figs. 

29a and 30a. We call this number s to denote a standard devia- 
tion ratio, and plot is as (JT = + s n_ particles/ cm 1 on Figs. 

P 

29a and 30a, for the total density and plasma density curves, 
s was calculated to be + 0.17. n_ j is used to denote any ordinate 
to the density plots and (T the standard deviation in particles/ 
cm^ . 

Now it is quite interesting to note that neither density 
corresponds to the Langmuir frequency within the confidence 
limits when we attempt to infer it from the total density plot 
However, when we explain our observed frequencies by correlating 
them with the plasma density curve, we note that these densities, 
n^ and ^ , fall within one standard deviation of this curve, if 
we consider that part of the plasma density curve at the point of 
disintegration of the beam density curve. 

The conclusions to be drawn from this analysis are brought 
out in the next section of this report. 

It is worthwhile to note that only the 472 Mcs oscillation had 
a large amplitude (2 ma of detector current) . The 520 Mcs oscilla- 
tion was much smaller in amplitude (0.01 to 0.02 milli-micro amp 
of detector current) . We associate the strong 472 Mcs oscillation 
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with the intermediate energy beam electrons (eV^ = 22 electron volts) 
and the 520 Mcs oscillation with the high energy beam electrons 
(eVb = 50 electron volts) . Note also that the density of this 
group of high energy beam electrons is about one order of magnitude 
smaller than either the plasma electrons or the intermediate energy 
beam electrons. 

Other high intensity electron oscillations were detected at 
higher currents much more readily than at 50 ma discharge current. 
Unfortunately, this was discovered so late as to preclude taking 
correlating probe data. The results of these measurements, for 
one probe position, are shown in tabular form below. 



Discharge Current 
(ma) 


Observed 

Frequency 

(Mcs) 


Approximate* 

Computed 

Frequency 

(Mcs) 


Intensity of 
Signal 
(Detector 
Current - amp) 


100 


703 


710 


3 x 10 -11 


105 


704.5 


720 


1.2 x 10~ 9 


110 


708.5 


730 


4 x 10 -7 


115 


71J 


737 


4 x 10” 7 


120 


713.3 


746 


4 x 10~ 7 


125 


717.5 


752 


6 x 10~ 7 



For this probe position no oscillation signals were detected 
below 100 ma discharge current. This suggests that a threshold 
effect for interaction may be present as a function of current 
density . 

* 

By interpolating in Fig. 14 using a pressure correction factor. 
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Fig. 40 is a plot of current density (J^) vs observed fre- 
2 

quency squared (f^) showing the linear relation of these values. 

The other type of oscillation found was a low frequency sig- 
nal of about 100 KcSo The ion frequency predicted by calculation 
from the density of the intermediate energy beam electrons is 
135 Kcs. Hence, these oscillations are believed to be ionic 
in nature. They existed throughout the length of the tube (for 
the 50 ma discharge at 6.9 microns pressure), maintaining 
essentially the same frequency; however, the intensity of the 
signal did vary with position in the tube. The results of these 
intensity variations are shown in Fig. 34, where I (intensity in 
volts) is plotted vs X. It is also significant to note that the 
oscilloscopic traces of these oscillations show a more fixed 
adherence to a steady fundamental frequency and waveform on the 
cathode side of the discharge where the high and intermediate 
energy electron beams, coming from the cathode, were in evidence. 

On the anode side of the discharge (in the region after which 
scattering has occurred) these oscillations show a wider spectrum 
of frequencies (varying between 90 and 110 Kcs) and the intensities 
of signals and waveshapes are qui Le variable when sampled in time. 
Fig. 41 is an oscillographic photograph of an oscillation pattern 
on the cathode side of the meniscus region. Fig. 42 is an oscillo- 
graphic photograph of an oscillation pattern obtained on the 
anode side of the meniscus region. 
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Visual and photographic observation of the discharge was 
made in effort to try to discover a meniscus, a scattering region 
and/or a diverging beam of electrons. No meniscus was observed 
at low currents; but, visual observation seemed to indicate a 
region of maximum disturbance in the beam, and at about the same 
location, a tendency for the beam to diverge from the shape dis- 
played in emerging from the cathode. The region in which this 
phenomenon appeared is indicated on Figs. 29 through 39 as the 
cross-hatched region. Figs. 43 through 45 are photographs and 
sketches of the discharge which show the beam, its tendency to 
diverge, and apparent meniscus formation. Note that meniscus 
formation appears to grow with increasing current. 

In attempting to correlate frequency measurements made at 
higher currents during the second data run (50 ma discharge at 
6.9 microns of pressure) with the current-density data of the 
first run (five different currents at 5 9 microns of pressure) , 
an interesting fact (overlooked as being trivial when it occurred) 
was noted. Referring to Fig. 14 for a discharge current of 125 ma, 
at a position ten mm from the cathode, a peak density occurs. The 
minimum of this gradient occurs at 12 mm from the cathode. Referring 
to the original probe characteristic curve corresponding to this 
point (Fig. 7, 125 ma discharge current), it was noted, that consid- 
erable irregularity of the ion portion of this curve was observed 
in a probe voltage range of -40 to -100 volts with respect to the 
anode. The phenomenon manifested itself in a step increase in 
probe current of about 0.004 ma which corresponded to a step 
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increase of about 0.4 volts. This distortion in probe character- 
istic lasted out to about - 100 volts at which point it tended to 
return to its normally expected shape, but with some noticeable 
change in slope. At the time of occurrence, the peculiarity was 
noticed and the curve retraced several times, each time giving a 
similar result, although the exact voltage at which the step in- 
crease in current occurred did not coincide exactly, it held to 
within a 20 volt region. Sweep speed of probe voltage was in- 
creased with no change in effect -- in fact the speed of sweep 
was increased to about 50 cps (vs the normal 0.05 cps) and the 
trace photographed on the oscilloscope. The same sort of 
phenomenon showed thereon except that this trace seemed to show 
oscillating signals riding on the probe characteristic from about 
-20 volts out to the maximum excursion of -130 volts. Further 
examination of the electron portion of the probe characteristic 
curve disclosed that this was the first probe characteristic curve 
(for this current) , in moving the probe from cathode to anode, in 
which a section having a constant slope in the positive current 
region (denoting a region separating two energy groups, ref. 

discussion of Waymouth’s work, pg. 23 , Section 2, Previous 

* 

Work) did not exist. 

A similar occurrence took place for the same discharge current 

at 10 mm from the cathode. This disturbance to the probe was less 

pronounced, amounting to a change in slope on the linear ip + vs V 

2 

curve and a distinct deviation from linearity on the ; vs V„ curve. 

P+ P 

See Fig. 46 which illustrates this. cf. Fig. 8, in which only 
one electron group is present. 



55 



Still considering this 125 ma discharge, it was noted at 
the time of taking the data that the position 10 mm from th< 
cathode was just barely on the cathode side of the diverging 
beam region (observed visually) of the discharge- The electron 
portion of this probe characteristic curve showed a constant 
slope region denoting two distinct groups of electron^ still 
present at this position in the tube. The position 12 mm from 
the cathode was inside, i.e., on the anode side, of the diverg- 
ing beam. The meniscus is believed to have been present about 
10 mm from the cathode. 

Similar occurrences were evident on the probe character- 
istic curves at 14 mm from the cathode for the 100 ma discharge and 
at 22 jiun (see Fig. S) from fhe cathode for the 75 ma and 125 ma, dis- 
charges at probe voltages of about -100 volts with respect to 
the anode. These points corresponded with a peak density 
gradient in each case (see Fig. 14) . 

In light of what was discovered on the second data run 
(50 ma discharge at 6.9 microns pressure) it now appears clear 
what happened in these cases during the first data run. These 
disturbances to the probe characteristic are associated with the 
same region in which (a) electron oscillations were detected, 

(see Fig. 29). (b) ion oscillation intensity reached a maximum 
and ion oscillation waveshape became most distorted, (see Tig. 34) „ 
(c) beam electron concentration largely disappeared and merged into 
one random distribution (see Figs. 23 and 24), and (d) either a 
maximum or a minimum density gradient was present (see Fig. 16) 
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The changes in slope as far negative as -100 volts ace 
believed to be caused by electrons acquiring enough energy in 
the interacting region in the vicinity of the meniscus (where 
oscillations are believed to originate) so thit they can overcome 
this amount of retarding field and become collected by the probe 
There are, it appears, distributions of electrons in these inter- 
acting regions having maximum energies as high a- LOO + electron 
volts. The step increase in current shown on Tig 7 is believed 
to be a monoenergetic group of electrons having an energy of 
about 50 electron volts with a variability of about + 10 electron 
volts when sampled in time. 

A distinction can be made with respect to the effect on the 
probe characteristic as it is affected by the density gradient, 

A large disturbance, i.e., a step increase in probe current s and 
a small disturbance, i.e., a change in slope^was noted where the 
position was associated with a minimum of the density gradient 
(corresponding with Emeleus T findings that oscillations reach a 
maximum intensity at the bottom of a density gradient) . Only a 
small disturbance, i.e., a change in slope of the probe character- 
istic curve, was noted where the position corresponded to a maxi- 
mum of a density gradient (in agreement with Emeleus* findings 
that oscillations are initiated at or near peaks in the density 
gradient) . 
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The effect of probe contamination was studied to determine 
whether cathode sputtering was causing the probt to ro'lect im- 
pure material and thus change its collecting properties. A 
curve was traced with a freshly cleaned probe J^by ohmic, heating 
to a cherry-red condition) then allowed to sit immediately in 
front of the cathode, perhaps 2 mm away, (the worse possible 
position), for a period of 15 minutes to one ha^f hour and then 
the same characteristic curve was retraced. A slight change in 
the characteristic was noticeable; however, when the interval of 
time was reduced to the same time as was used for making a set of 
data curves, the effect of probe contamination on the character- 
istic curves was entirely negligible. This is probably uttcibut 

able to the very low evaporation rate of material from the cathode, 

2 

given as less than 0.8 micrograms of barium per cm per hour by 
the manufacturer. 

Miscellaneous tests were conducted to determine the effects 
of electronic errors associated with the measurement of probe 
characteristic curves. Probably the best estimate of the relia 
bility of the probe characteristic curves was obtained bv allow 
ing the curve to be retraced 15 to 20 times to note the drift 
tendency due to the amplifiers and associated electronic circuitry 
Results of these measurements showed that the maximum deviation to 
be expected was of the order of + 0.002 ma of current for electron 
currents of the order of 0-1 rna. Curves were also run after an 

elapsed time of days and no significant change in characteristic 
was noticed except when the gas was obviously dirty because of 
long periods of operation (up to several weeks) . 
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60 Conclusions. 



The results of our second data run (50 ma discharge at 6 9 
microns pressure in argon) as shown in Fig. 29 tend to substanti- 
ate that plasma electron oscillations originate at the peak of a 
density gradient and grow in amplitude down a density gradient as 
originally proposed by Allis (30) and observed experimentally by 
Emeleus and co-workers (33) . 

Electron oscillations were detected only at positions in the 
tube immediately after divergence of the electron becim had occurred 
These oscillations could be associated with a peak density gradient 
and the observed frequency of these electron oscillations bore a re 
lationship to the predicted Langmuir frequency when correlated with 
the plasma density (vs the usual total density) . The reason for 
this will become clear later in discussing a possible mechanism 
for oscillation excitation. 

Ion oscillations were detected throughout the length of the 
tube with a fundamental frequency of about 100 Kcs. Fig. 34 shows 
that the intensity of these oscillations fell to a minimum just 
prior to the point where the total density gradient reached a 
maximum (see Figs, 29 and 34) and then the oscillation intensity 
rose rapidly to a maximum where the total density gradient Wcas a 
maximum. This region is presumed to be that where particle inter- 
action is greatest and where energy exchange is taking place and 
thus exciting electron oscillations. This region is outlined by 
cross-hatching on Figs. 29 Through 39. The fact that the photo- 
graphic traces of the ion oscillations across this region., in 
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going from cathode to anode, show considerable regularity on tho 
cathode side (see Fig. 41) and more intensity and more irregularity 
in wave shape and frequency on the anode side (see Fig. 42) 
indicates that in this region, strong m.icrofields are present 
which are disturbing the ion oscillations. The change in intensity 
of the ion oscillations indicates that ions are taking part in the 
energy exchange mechanism by losing energy as they move out of the 
turbulent region toward the cathode. 

Disturbances to the probe characteristic curves noted on the 
first data run being associated with maximum or minimum density 
gradients points strongly to the conclusion that strong fields 
and oscillations were present, and considerable energy exchange 
was taking place, at these positions in the tube. Disturbances 
out to -100 volts suggests that some electrons were acquiring 
better than 100 electron volts of energy These fast electrons 
must be generated in the interacting region since the maximum, 
energy that could be gained from the discharge tube would be 
of the order of the tube voltage drop, about 20 electron volts. 

Maximum disturbances were noted at positions corresponding 
to a minimum of the density gradient in agreement with Allis v 
contention that oscillations grow in intensity down a decreasing 
density gradient. Smaller disturbances to the probe character- 
istic curve were noted at positions corresponding to a maximum 
of the density gradient in accord with Allis ? idea that plasma 
oscillations originate at, or near, a peak density gradient. At 
higher discharge currents, multiple density gradients, (i.t., two 
peaks) were present (see Fig. 14) . Probe disturbances were 
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associated with each of these. This fact sugge^tb that a more 
complex interaction is taking place with perhaps the form it ion 
of two such interacting regions, one following the other. 

Examination of the distribution functions shown in Figb. in 
through 28 show that two groups of electrons are present on the 
cathode side of the interacting (meniscus) region. On the anode 
side of the interacting region, only one group of electrons is 
present. Further, by employing WaymoutHs analysis of the probe 
characteristic curve, a third higher energy group of electrons is 
also found to be present, but these electrons ar^ either much less 
numerous or of too high an energy to show on the energy distribu- 
tion function. We have earlier denoted these higher energy groups 
as an intermediate energy group of beam electrons and a high energy 
group of beam electrons, respectively. 

These beam electrons emerge from the region in front of the 
cathode and interact with the plasma electrons at a region in the 
tube where visual divergence of the beam and the aforementioned 
disturbances in the probe characteristic curves occur 

The evidence of both data runs has suggested to us that the 
mechanism exciting the oscillations is associated with these beam 
electrons, because all of these associated phenomena are indicative 
of a region of instability in which local fluctuating microfields 
of high strength would be likely to exist. 

In the light of this evidence, the authors believe that the 
beam electrons, upon reaching the plasma, are strongly interacting 
with the plasma (when beam current densities are adequate) 9 pro- 
ducing local instabilities which excite plasma oscillations. 
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Ion oscillations may well play a role in the excitation process , 
or may only be a by-product of the electron beam excitation 
mechanism. 

The theory of such a mechanism has been developed by Akhiezer 
and Fainberg and is discussed by Kharchenko, Fainberg, Nikolayev, 
et al., (31). The important points in their discussion are sum- 
marized in Section 2, Previous Work, page 26. 

As previously mentioned, there are two possible types of 
instabilities, convective and absolute. In our discharge, con- 
vective instability is believed to have occurred since oscillations 
were detected only in a very small region associated with the 
initiation of the instability, and none were detected on the anode 
side of the point of initial detection where one would ( xpe.et the 
oscillation intensity to increase as the oscillation moved down a 
density gradient. Possibly the spread in velocities as the beam 
electrons become randomized in their initial interaction with the 
plasma tends to remove the instability as suggested by Kharchenko, 
et al., (31). This may be the reason that oscillations are not 
detected other than in the region of initial scattering. 

Perhaps as a bit of speculation, we might mention that the 
possibility of modified "runaway" electrons (31, 1, 36) being 
generated and initiating the instability should not be entirely 
ignored. Although "runaway" electrons have not been observed 
(to our knowledge) in discharge tubes such as ours, conditions 
suitable for their production appear to be met in our tube and 
some of the phenomena observed might well be attributed to their 
presence. For example, the sudden termination of the instability 
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and damping out of oscillations are what one would expect if the 
excitation was produced as a result of "runaway" electrons. 

Or possibly, electrons traversing the interacting region may 
be given a pulse of energy by the oscillations and emerge as 
"runaway" electrons. These "runaways" would not interact with 
other particles because their collision prol V is 

less at their greater energy, thus, they would completely emerge 
from the plasma and finally give up some of their excess energy 
by collision with walls of the tube. 

One other point is worthy of discussion at thio timo On the 
second data run, oscillations were detected at higher discharge 
currents having signals of considerably greater amplitude than 
those observed at 50 ma discharge current. All efforts at detecting 
oscillations at discharge currents lower than 50 ma were unsucce.s 
ful. Examination of the density plot from the first data run 
(Fig. 14) shows that no sharp gradients exist for the 25 and 50 
ma discharges. Sharp gradients do exist for 75 ma and above. In 
the light of this observation and the experience obtained on the 
second data run, it is believed that current density (in order, to 
produce a sufficiently strong interacting beam) , and pressure are 
important in establishing the instability and. hence oscillations. 

In particular, had we operated at slightly higher pressures 
and/or slightly higher currents, it is believed tha T oscillations 
would have been found more readily and that: a closer approach to 
an ideal plasma (i.e., having a Maxwellian distribution of energi 
would have been obtained since a smaller mean free path would hav> 
produced more randomization of electron energies. 
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The possibility also exists that we were cperiting in a 
current density - pressure region which was just adequate to allow 
a beam to form and move some distance before interacting with th( 
plasma. Operation at higher current densities* has shown us th.it 
the meniscus region moves toward the cathode as urgent ii. increased 
Further experimentation should reveal the optimum current tensity 
pressure relationship to use in attempting to isolate the excitation 
mechanism. 

The results of this investigation are based upon only limited 
data, some of which is questionable in magnitude perhaps 10 to IS 
per cent in error in ordinates of the energy distribution function:) 
because of the equipment troubles at the end of our work, however, 
all data were qualitatively good, at least. Much more diti need 
to be collected before definite, positive conctusions may be drawn. 
Nevertheless, the correlation between all of the evidence pr* sen ted 
is quite convincing and seems to substantiate the tentative con- 
clusions which have been drawn. 



•j^ 

The area of the cathode used in this work is If cm £ ; h no 
dividing any of our discharge currents by if will giv> the current 
density . 
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7 . Recommendations . 



In keeping with our findings and observations, the following 
recommendations for further work are made: 

(a) Continue probe studies and oscillation detection ar 
higher pressures and/or higher currents, perhaps 
varying pressure as a parameter to determine changes 
in the form of distribution functions. 

(b) Establish the pressure-current density range in which 
oscillations are detected most readily, then attempt 
to find a simple oscillation-density gradient rela- 
tionship which will present the least complicated 
structure initially. 

(c) Once suitable oscillations are found, make probe 
studies and frequency measurements of both ion and 
electron oscillations to determine where in the tube 
they arise and attempt to correlate these measure- 
ments with visual observation, energy distribution 
functions, and perhaps also some spectroscopic 
observations . 

(d) Study the region in which oscillations appear to 
arise in great detail; i.e., make probe and oscilla- 
tion measurements around this region utilizing very 
small increments of movement of the probe, say O.S 
mm or less, in an effort to isolate the mechanism 
of excitation. 
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Figure 1 

General Experimental Arrangement 
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Figure 2 

I r- ‘’0 Characteristic Amplifying and Recording Eaui 
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Figure 3 

Power Supplies and Probe Circuitry 
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Figure 4a 

Discharge Tube 
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Figure 4b 
Discharge Tube 
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Figure 11. In i vs. V 
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Figure 12. In i vs. V 
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Figure 17 a. 

Probe Characteristic Curves 
i vs . V 
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Figure 17 b. 

Probe Characteristic Curves 
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Figure 18 a. 

Probe Characteristic Curves 
i vs . V 
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Figure 29. Electron density as a function of Cathode ° Fiobe position. 
Obtained from energy distribution functions. 

• -- Total density 

O -- Density of plasma (low energy group) 

X -- Density of beam electrons (intermediate energy groip) 





Figure 30. Total electron density 
position. Obtained from the slope 



as unction of Cathode - Probe 

of i vs. V curves. 
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Figure 31, Relative total electron density as a function of Cathode - 
Probe position. Obtained from In i vs. V curves ne&~ space potential. 




Figure 32. Relative electron density of beam electrons as a function 
of Cathode - Probe distance. 

0 -- Intermediate energy group 
X -- High energy group 
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Figure 33. Average energy of electrons in electron volts as a function 
of Cathode - Probe distance. Obtained from energy distribution functions* 
Intermediate energy group 
— Low energy group 




Figure 34. Ion oscillation intensity as a function of Cathode - Probe 
distance . 
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Figure 35. Electron temperature in degrees Kelvin as a function of 

Cathode - Probe distance. Obtained from In i vs. V curves. 

P P 

X -- High energy group ) In i vs, V curves corrected for ion current. 

• -- Low energy group } and Fesolved Snto a high and low energy group , 

0 -- Low energy group ) Uncorrected linear portions of i vs , V curves 

used for these data. ^ 




170 



o 



150 







130 



Figure 36. Average electron temperature in degrees Kelvin as a function 
of Cathode - Probe distance. Obtained from energy distribution functions. 

X -- High energy group, T 

e H 

O -- Low energy group, T g 



95 



Argon 




C A 

X (mm) 

Figure 37. Space potential in volts as a function of. Cathode - Probe 

distance. Obtained from maximum of first derivative of i vs. V 

P P 



Argon 




X (mm) 

Figure 38. Space potential in volts as a function" of Cathode - Probe 
distance. Obtained from the first deviation from linearity of In i 

vs. V curves. 
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Figure 39. Electric field strength as a function of Cathode - Probe 
distance. Obtained from data of Figure 37a 



CM 

^ S 

u o 



O) 

G 

<U 

Q 



C 

<D 

U 

U 

3 

U 



a 

§ 



CM 

I 



T3 




Figure 40. Current density as a function of the square of the observed 
frequency of electron oscillation for the probe located immediately 
after the region of suspected maximum interaction. 



97 



Fig 41a. Several oscilloscope traces 
of ion oscillations. 




Fig 41b. Oscilloscope trace of ion 
oscillation when sampled in time. 




Fig 41c. Oscilloscope trace of ion 
oscillation when sampled in time. 



Fig 41. Oscilloscope traces of ion oscillations on the 
cathode side of the meniscus region. Discharge current 
50 ma. Pressure 6.9 microns of mercury. Cathode-anode 
spacing 40 ran. Probe position 12 mm from cathode - 
oscilloscope sweep speed 5yi/sec/cm. Amplitude 1 volt/cm. 
Mote cm ller amplitude, more fixed waveform and more steady 
frequency compared with Fig 42. 
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Fig 42a. Several oscilloscope traces 
of ion oscillations. 




Fig 42b. Oscilloscope trace when 
sampled in time. 




Fig 42c. Oscilloscope trace when 
sampled in time. 



Fig 42. Oscilloscope traces of ion oscillation* 
on the anode side of meniscus region. Cathode-anode 
spacing 40 mm. Probe position 20 mm from cathode. 
Discharge current 50 ma. Pressure 6.9 microns of 
mercury. Oscilloscope sweep speed 5 ^sec/cm. 
Amplitude 1 volt/cm. Note the greater amplitude, 
distortion of wave form and spreading of fundamental 
frequency compared with Fig 41. 
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Figure 43 a. Photograph of 50 ma discharge 
in argon at 6.9 microns of mercury. 




X ( mm ) 

Figure 43 b. Sketch of 50 ma discharge in 
argon at 6.9 microns of mercury. 



Figure 43 
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Figure 44 a. Photograph of 120 ma discharge 
in argon at 6.9 microns of mercury 
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Figure 44 b. Sketch 
argon at 6.9 microns 



of 120 ma discharge 
of mercury. 



in 



Figure 44 
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Figure 45 a. Photograph of greater than 150 ma 
discharge in argon at 6.9 microns of mercury. 
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Figure 45 b. Sketch of greater than 150 ma 
discharge in argon at 6.9 microns of mercury. 



Figure 45 
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and we use the trapezoidal approximation to the area, then 
r =1 6 

N 1 i A eV £N0(eV) 2 + N0(eV) t 



oC«*2 



N 



N x i 1.722 x 10 



13 



(AeV) 



liv^i ' 1 

J_ 2 eV»2 



to 

[ oC«28 

ZH N0(eV) oC + iN0(eV) 

o(«l6 



'* c] ’ 



where i 



it 



d 2 i. 
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+ y * v s 

ct« 

Since we are usj.ng l/3 i" , we multiply by 3 and factor '“'P 

out 10 ”° from V 2 l" terms and our computational formula becomes: 

<< =16 



dv: 



N x i 1.722 (3) (2 eV) x 10 



L y2/2 \5 



y cC 




where y 



^i" x 
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It is convenient to express this computation in tabular form 
in which a running integral may be set up and the integral 
eV 2 

J V N0(eV) d(eV) computed at the same time for the purpose 
eV l 

of evaluating the average energy, eV, of a group of electrons. 



Figure 47. Example of electron density computation and average 
energy computation from energy distribution function. 
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f N0(eV) d (eV) i" d(eV) 



By the same process, N£ and eV may be computed. 

10 H 

Then, N " N l + N 2 m j N 0< eV ) <i(©V) 



These computations are shown in a tabulation on the 
next page in Figure l*7b. 



Figure k7 . Example of electron density computation and average 
energy computation from energy distribution function. 
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